Packing defects and the width of biopolymer bundles 
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The formation of bundles composed of actin filaments and cross-linking proteins is an essential 
process in the maintenance of the cells' cytoskeleton. It has also been recreated by in-vitro experi- 
ments, where actin networks are routinely produced to mimic and study the cellular structures. It 
has long been observed that these bundles seem to have a well defined width distribution, which has 
not been adequately described theoretically. We propose here that packing defects of the filaments, 
quenched and random, contribute an effective repulsion that counters the cross-linking adhesion 
energy and leads to a well defined bundle width. This is a two-dimensional strain-field version of 
the classic Rayleigh instability of charged droplets. 

PACS numbers: 



Filamentous biopolymers, such as F-actin, have the 
ability to cross-link into a variety of bundles and net- 
works, forming the cytoskeleton. The distribution of the 
radii of cross-linked actin bundles is a basic characteris- 
tic that determines the mechanical properties of the cy- 
toskeleton. The thickness of bundles similarly determines 
the mechanical properties of artificial networks that form 
in-vitro [l], Q ■ Recent in-vitro experiments 11 1 indicate 
a broad distribution of radii, with a distinct peak, while 
simple equilibrium theory would predict a global phase 
separation and formation of a single bundle of infinite 
width and length (in an infinite system) . Several pos- 
sibilities have been proposed to explain the observed dis- 
tribution; When the bundles form due to multi-valent 
ions, electrostatic interactions are unbalanced and can 
lead to a selection of an equilibrium finite radius 0]. A 
recent study attributes the finite size of the actin bundles 
to the inherent chirality of these filaments [7| , and indeed 
this may be the dominant effect when the filaments are 
linked by small multi-valent ions [1, Q. We will deal 
here with neutral systems that are strongly chemically 
cross-linked by larger linker proteins 0, U^|Tl| . In these 
cases the chirality of the individual filaments may play 
a lesser role, and we will treat the filaments as achiral. 
Our model proposes therefore a complimentary mecha- 
nism for finite bundle widths, which applies to achiral 
bundles. We present here a model where quenched dis- 
order, in the form of twist-defects (Fig.l), leads to the 
selection of a quasi-cquilibrium finite radius. We call this 
a quasi-equilibrium since the defects in the bundle can be 
annealed away in principle, but this is highly unlikely for 
long filaments and strong cross-linking, and the system 
is therefore dynamically arrested in a meta-stable con- 
figuration. There is recent experimental evidence that 
in-vitro actin bundles do indeed have such twist defects 
3 . We wish to consider the implication of such defects 
for the width distribution of in-vitro actin bundles. We 
will consider these knots to be in a form of quenched, i.e. 
static, and random disorder. 

Consider a bundle of cross-linked filaments, as shown 
in Fig.l. The filaments are assumed to be tightly cross- 



linked to all their neighbors, so that they form a solid 
structure. If all the filaments are perfectly aligned then 
they form a hexagonal close-packed crystal. But in the 
process of aggregation twist-defects can form, where two 
filaments attach to the growing bundle in a different rela- 
tive orientation along their lengths (Fig. la). If these two 
filaments stay on the surface of the bundle and have a 
single such twist along their length, then this twist can 
slide to the ends of the filaments and relax away. When 
the filaments are long enough it is likely that they each 
form such twist-defects with other neighboring filaments 
in other places along their lengths, and the whole bundle 
is in fact knotted (Fig. la). The knots can not be easily 
relaxed now by sliding them to the end of the filaments, as 
they are all entangled. The meandering of the filaments 
in the bundle due to the twist-defects may be described as 
a random walk of the filament, where the average number 
of defects along its length is: Ndef ~ {L/Lp)Pdef, where 
L is the overall length of the filament, and each twist- 
defect extends over a length Lp which is of the order of 
the persistence length of actin in the bundle (Lp ~ 16/im 
for free actin filaments [l^l), and Pdef is the probability 
of exciting a defect. 

When actin bundles attract each other, through the 
strong adhesion of cross-linking proteins, they begin to 
aggregate laterally. We will consider here the case of 
strong cross-linking and long filaments, so that the inter- 
action term dominates over the entropy term in the free 
energy, which we will neglect from now on . The adhe- 
sion energy gain is driving the aggregation, and since the 
filaments on the surface have less adhesion energy (less 
neighbors), there is an energetic drive to increase the 
bundle thickness. The adhesion energy per unit length 
can be simply written as 



E^nd « -TT{R/a)ebiiR/a) - 2) 



(1) 



where is the adhesion energy per unit length, between 
the filaments due to the presence of the cross-linking pro- 
teins. The first term represents the bulk adhesion energy 
(negative) and the second term represents the surface 
energy (positive). This energy functional gives a critical 
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radius oi R = a {a is the radius of the individual fila- 
ments and surrounding cross-linkers, which is typically 
0, [Hi ~ lOnm). beyond which the bundle grows to infi- 
nite width. The observed widths of actin bundles seem 
to be rather well-defined and despite the mobility of 
the filaments the bundles only grow by longitudinal ag- 
gregation and do not continue to thicken through lateral 
aggregation. This observation is therefore at odds with 
the equilibrium theory described by Eq. ([!]), and is our 
motivation for looking at the effect of packing defects. 

We now treat the case of twist-defects inside the bun- 
dle. Each defect involves an increase in the local energy 
due to several terms (Fig.l): (i) Loss of adhesion energy 
due to broken cross- linking bonds at the site of the twist, 
(ii) the elastic energy of the twisted actin filaments, (iii) 
the elastic energy of the deformed hexagonal lattice of fil- 
aments around the defect in the twist plane (Fig.l). Out 
of these energies the first and third can be different if the 
defect is inside the bulk of the bundle or close to the sur- 
face. Furthermore, it is clear that these two terms arc in 
fact smaller at the surface, so the energetic cost of a fil- 
ament is higher inside the bulk of the bundle. A uniform 
distribution of static defects will therefore cost a higher 
energy the thicker is the bundle. We can already see how 
the defects may compete with the adhesion energy and 
drive the bundles to a finite preferred width. 

Let us treat the bundle as being approximately uniform 
along its length, so that we can treat the defects in a two- 
dimensional circular cross-section. The strain field in the 
surrounding bundle due to a localized defect depends on 
the order of the defect, i.e. a monopole or a higher multi- 
pole, and is maximal in the maximal twist-plane (Fig.l). 
A twist of two filaments increases the local volume of the 
filament packing (Fig.l), so has a monopole component 
(dilation). At the same time the twist involves also a 
quadropole component, as shown in Fig.lc. The single 
isolated monopole defect gives rise to a strain field that 
decays as u ~ 1/r away from the defect, and the result- 
ing strain energy has a logarithmic divergence with the 
radius of the bundle. The contribution from such a single 
defect is therefore negligible compared to the term in 
Eq.©, and the equilibrium stays at an infinite width. 

We are therefore led to consider a uniform distribution 
of many defects. Before dealing with the complications 
of the long-range strain field of the monopole, let us treat 
the case of a highly localized strain-field, such as for the 
quadropole. In this case we can treat the strain energy 
in the bundle around the defect as part of the defect core 
energy Ec, that includes the broken cross- linker bonds 
and twist of the actin filaments. 

A uniform distribution of defects, of density p = 1/i^, 
will simply add Edefect = pirR^Ec to the energy of 
Eq.©, without any qualitative change, i.e. the system 
still has minimum energy for an infinitely wide bundle. 
Let us assume that close to the rim of the bundle defects 
are not strongly trapped, and can relax. Such a process 
may occur due to higher mobility of the filaments on the 
bundle surface, which leads to effective " surface- melting" 



of several layers of filaments, thereby annealing any de- 
fects. The thickness of this annealed layer is denoted by 
A (Fig. 2a). The energy of the defects (per unit length) 
in this system is therefore given by 

Edefect ~ - >^fEc (2) 

The total energy per unit length is now Etotai = 
Ebind + Edefect ('Eas lll2[) . Wc find that there is a crit- 
ical value of the defect energy E* = eb(L/a)^, above 
which the system has an equilibrium configuration of fi- 
nite width, given by 

Ro^sheii = {EeX - aE:)/{E, ~ El) (3) 

Note that as the defect core energy increases above its 
critical value, the bundle radius shrinks, for a fixed 
density. There is also a critical value of the width of 
the annealed layer, given by: Ac = aE*/Ec, at which 
the equilibrium bundle size shrinks to zero. Above E* 
the system is dominated by the energy of the defects, 
while below it the binding energy dominates. In the 
limit of Ec ^ eb,E*, we find that the equilibrium ra- 
dius approaches the thickness of the defect-free layer 
Ro, shell — > A, while in the limit Ec — > E* the equilib- 
rium radius diverges. 

When the width of the defect-free outer shell is given 
by the typical inter-defect distance, i.e. A = L, we get an 
equilibrium configuration of finite width, given by (EqjS]): 
Ro,sheii = L{Ec - ebL/a)/{Ec - E*), and a critical inter- 
defect separation of: Lc = a{Ec/eb)- 

We now return to the problem of a uniform distribu- 
tion of strain-monopoles. These can be treated as a two- 
dimensional gas of charges, all of the same sign, interact- 
ing via logarithmic repulsion [T^] (Fig. 2b, bottom). Note 
that since the twist-defects cause a local dilation, they 
all have the sam sign and repel each other. The prob- 
lem therefore resembles a two-dimensional version of the 
famous Rayleigh- instability of charged droplets [l^ . In 
our case the twist-defects behave as elastic charges and 
their mutual repulsion breaks the infinitely thick bundle 
into bundles of finite radius. 

We can estimate the overall strain energy (per unit 
length) of the system of uniform monopoles by calculat- 
ing the interaction of a central charge with a uniform 
distribution of surrounding charges, assuming that its 
nearest neighbor is a distance L away (Fig. 2b) 

Emon ~ k\n{L/a} + 2kp J \n{r/a)(fr (4) 

= A: In (L/a) +2k^ - ^ - In {L/a) + ^ In {R/a)^ 

where k is the elastic modulus of the bundle which is of 
order Sb- The first term in Eq.© is the energy to create 
the defect at a distance L from its nearest neighbor, and 
the second is the interaction energy with a uniform distri- 
bution of surrounding defects. We find that the dominant 
term in R is of order {R/L)^\n{R/a). Since we do not 
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attempt an exact solution of the complex strain-field in- 
side a cylindrical bundle with defects, we will continue 
with a scaling analysis using this term, i.e. we will ap- 
proximate the strain-field energy per unit length due to 
the defects as 



The width of the distribution can be estimated from 
the above calculation, using 



(5) 



where k' combines the effective stiffness of the bundle 
and various geometric factors. This energy dominates 
at large R over the binding energy (Eq[T]), due to the 
In (R/a) factor. We therefore have a situation now where 
the bundles always have a finite equilibrium radius. 

The total energy per unit length is now Etotai = 
Ebind + Edefects (Eas lll5| ). and has a global minimum 
at the equilibrium radius, given by 



Ro.mono/a = exp ((l/2a^/9a) - 1/2) 



(6) 



where a = k' jsh- We find that as the density of defect 
increases the radius decreases, while it also decreases with 
the stiffness of the bundle a. Note that the equilibrium 
radius has an exponential dependence on the density of 
defects, and is therefore predicted to be very sensitive 
to this parameter. The calculation we have given above 
are all appropriate for a system with an infinite reservoir 
of actin filaments and cross- linking proteins. In a closed 
system with a finite number of available filaments and 
cross- linkers we have to conserve the overall area of actin 
bundles, which should not change the overall behavior we 
described. 

If the defects are interacting strongly with each-other 
through their quadropole component of the strain, then 
they can arrange in string- like aggregates (Fig.le), in the 
maximal twist-plane. These aggregates can form if the 
density of defects is high and they are mobile enough 
to move inside the bundle and aggregate. Alternatively 
such aggregates of defects may form as the filaments join 
the bundle at the surface, since a defect is more likely to 
form close to an existing defect; at the right orientation it 
has a lower strain energy (Fig. Id). Such aggregations of 
defects will modify the strain-field and the contribution 
of the defect energy, but is beyond the scope of this paper. 

Let us now compare our calculation with recent in- 
vitro measurements of the distribution of actin bun- 
dle radii The size distribution of bundles 



can be approximated to follow from a simple Boltz- 
man distribution of the bundle energy: PiR) oc 
cxp{—{Etotai{R) — Eo)/kBT. Since the bundle energy 
has a quadratic minimum at Rq (Eas l3l6| ). this gives a 
Gaussian distribution of radii centered around the equi- 
librium radius. The distribution of the number of actin 
filaments inside the bundles N is simply the distribution 
of bundle areas, N cx [R/af', which we we predict to 
have roughly an exponential behavior at large N since 
both Eas. (|2l5| ) are quadratic with R (for the monopoles 
there is a slight logarithmic correction). This is in good 
agreement with the observed distribution (Fig. 3). 



(Ai?2) 



knT 



\Rq 



ksT ksT 
2p{Ee-E*)' 'Wp 
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where the first is for the empty shell case fEas l2l3p . and 
the second is for the uniform field of monopoles fEqs l5l6p . 
We find that as the density of defects increases the width 
of the bundle sizes distribution decreases, and the distri- 
bution becomes tighter around i?o- 

In Fig. 3 we plot the calculated distribution compared 
to the experimental measurement at a ratio of 1:5 
Fascin cross-linkers to actin in the solution. In these 
fits we have to choose the unknown binding energy, and 
then fit the value of the parameters so that we get the 
position of the peak of the distribution and its width 
(Eql?]). There is some freedom in the choice of these 
parameters; for the empty-shell model we choose the val- 
ues; £{, = lOfcsT/a which is of the order of the known 
adhesion energy per typical actin cross- linker and 
a? p = 0.04 which corresponds to a rather dense array of 
defects and is quite arbitrarily. By fitting to the experi- 
mental distribution, we then fix the values of A/a = 1.007 
and Ee = 250fcsT/a. 

For the uniform-monopoles the width of the observed 
distribution fixes the value of pk' , and we take p = 
0.1 and k' = 0.03kBT/a, which then fixes the value of 
Eb — O.OS/csT/a, by fitting to the peak location. We find 
that for a good fit we need a very small value for the 
binding energy per unit length. The reason for this may 
arise from the fact that our two-dimensional monopoles 
represent twist-defects which extends over a length of 
the order of the persistence length of actin Lp [l^ . This 
may be much longer than the effective thickness of the 
maximal twist-plane Lmono (Fig.l). If Lmono ~ a then all 
the energy terms are re-scaled by a factor of Lmono/ Lp ~ 
10-2 - 10-3. 

It has been further observed that the average thickness 
of the bundles increases with the concentration of cross- 
linking proteins 0, [ll|. Within our model the radius 
of the peak of the distribution, at Ro^sheii or Rq, 
depends only on the density of defects p if both the bind- 
ing energy ei, and the defect energies Ec,k' depend in 
the same way on the concentration of cross-linking pro- 
teins (see Eqs l2l3p . We can therefore propose that as 
the concentration of cross-linking proteins increases the 
density of defects decrease, leading to an increase in the 
average radius of the bundles; both the peak position 
and the width of the distribution (Eq[7]) increase with 
decreasing p. The more weakly bound bundles appear 
therefore to aggregate in a "messier" fashion, resulting 
in more twist-defects. The probability of exciting a de- 
fect Pdef ^ exp {—/S.E /ksT) per unit length of filament, 
is higher when the energy barrier to create a defect at 
the bundle surface Ai?, decreases with the decrease in 
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the concentration of cross linkers. This prediction may 
be checked in future experiments. 

In a Uving cell the actin filaments form bundles in a 
variety of forms. Some bundles are inside the bulk of the 
cell cytoplasm, and their formation may therefore resem- 
ble that of in- vitro bundles, and consequently their width 
may be limited by the appearance of defects. Another 
form of actin bundles appear in the core of stereocilia, 
and seem to be packed in a perfectly regular, hexagonal 
lattice These bundles do not form by the lateral 

aggregation of preexisting long filaments as in the bulk 
in- vitro experiments; all the filaments in such a bundle 
polymerize together and elongate in synchrony from the 
stereocilia tip. The polymerization of the filaments at 
their ends is promoted by tip-complex proteins. If the 
relative locations of the nuclcation sites within this tip 
complex are maintained in a solid-like static structure, 
then the filaments which elongate from these nuclcation 
sites will form a perfectly packed bundle. The width of 
these bundles is therefore not controlled by the random 
process of defect formation, but is very tightly controlled 
by the cell [H. Indeed such bundles contain thousands 
of actin filaments, much thicker than the in- vitro bundles 



(Fig.3). 

We conclude that the aggregation of long filaments by 
strong cross-linkers is likely to produce packing defects. 
Such quenched and random defects can lead to an effec- 
tive pressure that limits the growth of the width of bun- 
dles, and may explain the observed peaked distribution 
with an exponential tail. It is a two-dimensional elas- 
tic version of the Raylcigh-instability of charged droplets 
[l^ . The rudimentary treatment given here should be 
improved in the future by a more rigorous treatment of 
the full three-dimensional problem of the complex strain- 
fields of a population of twist-defects inside a bundle that 
is free to twist and bend as a whole. 
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FIG. 1: (a) Schematic picture of a single twist-defect (red) 
inside a bundle of radius R and overall length L. The twist 
extends over a length of the order of the persistence length of 
actin Lp. Several twist-defects along the length of the bundle 
cause the filaments to form a knot. Far below or above the 
maximal twist plane (horizontal dashed line), the filaments 
are arranged in a two-dimensional hexagonal array (b). (c) 
In the maximal twist-plane the defect (red) creates a strain 
field around it (arrows) , and pushes the surrounding filaments 
from their perfectly hexagonal array (dashed circles) . (d) The 
quadropole strain field around a defect lowers the energy to 
create a near-by defect with orthogonal orientation (arrow), 
and can result in a string- like aggregate of defects (e). 




FIG. 2; (a) Schematic picture of the empty-shell model; 
uniform population of defects (red pairs) at areal density 
p = and core energy per unit length Ec, with an empty 

outer shell of thickness A. (b) Schematic calculation of the 
uniform field of monopoles; a single monopole (red pair in the 
center) interacts with a uniform distribution of neighboring 
monopoles of density p (red ring), beyond a minimal sepa- 
ration radius L. This system is the elastic equivalent to a 
two-dimensional electrostatic repulsion problem (bottom). 




FIG. 3: Experimental distribution of the radius of actin fila- 
ments inside bundles (circles) [ill ], for Fascin/ Actin ratio of 
1:5. The calculated distribution is given by the solid line for 
the empty-shell model (Eq[2ll and dashed line for the uniform 
monopoles (Eq[SJ . The inset gives a semi-log plot of the dis- 
tribution. 



